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Planetary nebula progenitors that swallow binary systems 

Noam Soker^ 

ABSTRACT 

I propose that some irregular ‘messy’ planetary nebulae owe their morphologies 
to triple-stellar evolution where tight binary systems evolve inside and/or on the out¬ 
skirts the envelope of asymptotie giant braneh (AGB) stars. In some eases the tight 
binary system ean survive, in other it is destroyed. The tight binary system might 
breakup with one star leaving the system. In an alternative evolution, one of the stars 
of the brook-up tight binary system falls toward the AGB envelope with low specifie 
angular momentum, and drowns in the envelope. In a different type of destruetion pro- 
eess the drag inside the AGB envelope causes the tight binary system to merge. This 
releases gravitational energy within the AGB envelope, leading to a very asymmetri¬ 
cal envelope ejection, with an irregular and ‘messy’ planetary nebula as a descendant. 
The evolution of the triple-stellar system can be in a full common envelope evolution 
(CEE) or in a grazing envelope evolution (GEE). Both before and after destruction (if 
take place) the system might lunch pairs of opposite jets. One pronounced signature 
of triple-stellar evolution might be a large departure from axisymmetrical morphology 
of the descendant planetary nebula. I estimate that about one in eight non-spherical 
PNe is shaped by one of these triple-stellar evolutionary routes. 

Subject headings: binaries: close — planetary nebulae 


1. INTRODUCTION 


Planetary nebulae (PNe) are believed by many to acquire their non-spherical morphologies 
through the interaction of their pr ogenitor with stellar or sub-stellar companions (e.g, see sum- 


erties of elliptical and bipolar PNe ( 

Sokei 

1998 

), and single stars cannot maintain a fast rota- 

tion over the evolution along the AGB (e.g.. 

Soker & Hamaz 1992: 

Nordhaus & Blackman 

20061; 

Garcia-Segura et al. 2014 

). Eollowing earlier studies, fe.g.. Bond & Eivio 

19901; 

Bondi 

2000l), 
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the b inar y shaping paradig m has gain a critical support in recent years (see reviews by Izijlstra 


liptict 

il and bipolar ] 

^Ne s 

ind a careful analysis of their morphologies, e.g. 

. Akras et al. (201.51). 

Aller et al. (2015a 

bl 

. BoffinI (2015). Corradi et al. (2015 

). De Marco et al. ( 

2015|), 

Douchin et al. 

(2015 

), 

Pang et al. 

c 

>015^ 

Hillwig et a' 

.1 (I2OI5). Ijones 

2015) 

Jones et al. 

(2015) 

Manick et al. 

(2015 

), 

Martinez Gonzalez et al. (2015). 

Miszalski et al. 

2015) 

Mocnik et a 

.(2015 

). Montez et al. 

(2015 

), limiting the list to papers from 2015. These close binary systems went through a com- 


mon envelope evolution (CEE). Some other non-spherieal nebulae around AGB and post-AGB 
stars might have been shaped by wider binary eompanions that did not go through a CEE (e.g., 
Eagadee et al. 2011 : Bujarrabal et al. 2013 : Van Winekel et al. 2014 : Deein et al. 2015 ). 


Despite mueh progress (e.g., Ilvanova et al.ll2013h . the CEE of these systems is not well un¬ 


derstood. Numerieal simulations cannot eject the envelope in a eonsistent and persistent manner. 


and spiraling-in ceases too early (e.g.. 

Sandauist et al. 

1998; 

Eombardi et al. 

2006 

De Marco et al. 

2011 

Passv et aklboil 

2012 

Ricker & Taamll2012 

). In light of these difficulties (e.g.. 

Sokeill2013l) 


it has been suggested that in many eases jets launehed by the more eompa ct eompanion, inelu ding 
main sequenee seeondary stars, faeilitate envelope ejeetion in a CEE te.g. . 180^12013 . 2014 ). 


If the jets launehed by the seeondary star remove the entire envelope outside the seeondary 
orbit as it spirals in, the formation of a CE is prevented. Instead, a grazing envelope evolution 
(GEE) takes place, where the binary system might be considered to evolve in a state of just entering 
a CE phase (ISokeiil2015l) . The companion grazes the envelope of the giant star, and both the orbital 
separation and the giant radius shrink simultaneously. Tidal interaction that leads to orbit shrinkage 
makes the GEE an alternative to the CEE in cases where jets efficiently remove the envelope. 
When the companion is massive enough to bring the giant envel ope to synchronizatio n, the orbit 
during the GEE does not shrink much, and might even increase. iGorlova et al. (I2OI2I) found that 


the companion to the post-AGB star BD+46°442 launches jets. The orbital period of this binary 
system is 140.77 davs. iGorlova et al.l(l2015b reported the detection of a collimated outflow from the 
companion of IRAS 19I35-I-3937, a post-AGB binary system with an orbital period of 127 days. 
Both these systems are highly compatible with the expectation from the GEE when the companion 
manages to maintain the AGB star in synchronization with the orbital period. I suggest that these 
systems evolved indeed through a synchronized GEE. 

The multitude of processes that are involved in the CEE and GEE phases, such as mass lose 
through the L 2 Eagrangian point and jets launching with an without precession, and the large 
parameter sp ace of binary system properties, ensure that each PN is ‘unique’ in its morphology 
(iSokenbOOTI) . As well, the multitude of processes place PNe on the crossroad of many astrophysi- 
cal objects, e.g., novae, symbiotic binaries, and massive binary systems such as rj Carinae, and of 

















































































































































-3- 


many generic astr ophysical eyolutionary phases, e.g., CEE, jet-inflated bubble s, and stellar merger 


sebrenko & Sokei 

2012[) . These pro- 


(IDe Mareoll2015h . PNe are directly connected also to type l a supernovae ( 

2 OI 5 I 1 andto intermediate luminosity optical transients (lEOTs; Soker & Kashi 
cesses and objects inyolve strongly interacting binary systems. 

On top of the large parameter space of binary system properties, a tertiary star can be added. 
In preyious papers mainly wide tertiary stars were mentioned. Namely, the secondary and tertiary 
stars orbit the AGB progenitor of the PN, and not each other. A wide stellar companion, to a central 
single star or a central binary system, at orbital separations of 03 ~ 10 — several x 10^ AU can 


cause t he PN to pos s es dep arture from axisymmetry and to haye an equatorial spiral pattern dSokei 
1994 ). Soker et al. ( 1992 ) proposed that the departure of the PN NGC 3242 from axisymmetry 
was caused by a tertiary star of a mass of M 3 IMq at an orbital separation of 03 ~ 4000 AU. 
In addition, a yery wide tertiary star, 03 > 10^ AU might form a small bubble inside the nebula 
(ISokeilll996l) . 


Engulfed binary systems were discussed in the past as well. lExter et al.l (120lOh proposed that 
the AGB stellar proge nitor of the PN SuWt 2 engulfed a tight binary system of A stars. They (see 
also lBond et al.ll2002l) proposed that PNe with a high density ring in the equatorial plane can result 
from triple-stellar eyolution. In that scenario the tight binary is engulfed by the primary star, and 
surviyes most or all of the CEE. It might merge later on as the stars of the binary system eyolye 
(Bond et al. 2002). 


In the present paper I add more to this rich variety of binary and triple evolutionary tracks by 
considering new evolutionary routes were instead of a single secondary star, a tight binary system 
is swallowed by the AGB or red giant branch (RGB) star. The typical initial orbital period of the 
tight binary system is Pzams,23 ~day—month, while the orbital period of the triple system (the 
tight binary and the primary star aro und their mutual cente r of mass) is Pzams,123 ~ 1 — 10 yr. 
Such systems are known to exist (e.g., Tokovinin et al.ll2006 ). This preliminary study sets down the 
foundations for further studies, including three-dimensional hydrodynamical simulations of triple- 
CEE and triple-GEE. In section |2] I describe the triple-stellar systems that are considered here. In 
section |3] I describe the different possible fates of such triple systems, and the general departure 
from symmetry of the descendant PNe. I summarize in section HI 


2. SETTING THE PROBLEM 
2.1. Expectation 

Eet us estimate the fraction of AGB stellar progenitors of PNe that interact with a tight binary 
system, rather than with a single secondary star. Consider main sequence stars that are expect to 
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evolve along the upper AGB, henee having a zero age main sequence (ZAMS) mass of Mzams,i — 
1 — SMq. Consider only those primary stars that have a main sequence stellar companion with a 
mass of M 2 < Mzams,i, with an initial semi-major axis azAMS,i 23 - The question is the following. 
What fraction of these systems have a tertiary star with a mass of M 3 < M 2 that orbit the secondary 
star in a tight orbit, azAMS ,23 < azAMS,i 23 ? I call this fraction 771 ^ 23 - 


In their review article lOuchene & KrausI (120131) give for stars with ZAMS mass in the range 
Mzams.i = 1-5 — 5Mq a multiplicity frequency of MF> 50% and a companion frequency of 
CF= 100 ± 10% (their Table 1 and figure 1) . For solar type stars the fraction of systems with 
n < 6 stars goes approximately as N{n) ^ fc2.5“", and Ri 25% of systems of solar type stars have 
n > 3 ( Duchene & KrausIboiB . section 3.1.5). For the mass range 1.5 — 5Mq that is more relevant 
for Galactic PNe, I take the fraction of systems with n stars to go as N(n) = kx~^ as well. I limit 
the systems for up to 6 stars. I find that to account for a MF> 50% and a companion frequency 
of CF= 100 the value is a: ~ 1.9 (taking x = 2 that gives CF= 90 does not change the results 
much). The fractions of systems with n stars, N{n) oc 1.9“", are as follows. Single star systems 
A^(l) = 0.49, binary systems N{2) = 0.25, triple systems, N{3) = 0.13, and then A^(4) = 0.07, 
N{5) = 0.04, and A^(6) = 0.02. These values give MF= 51% and CF= 


Let us first consider only binary and triple systems. In the triple-stellar systems the tertiary 
star can orbit the secondary star, the primary star, or both. We now consider the first possibility 
which can lead to a case where an AGB star swallows a tight binary system of two main sequence 
stars. I crudely assume that in about third of all triple-stellar systems the initial more massive star is 
at a larger orbital separation from the two lighter ones. This is a fraction of ~ N{3)/3 = 0.043 of 
all systems. The rest, ~ 0.09 of all systems, are triple-stellar systems that are not the triple systems 
consider now. The fraction of the systems considered now to all cases with binary interaction under 
these assumptions is 


tight (M 2 , M 3 ) sytems ^ N{3)/3 ^ 0.043 

all multiple systems N{2) -f N{3) 0.25 -f 0.13 


Three considerations increase the fraction of PNe shaped by triple-stellar systems relative to 
the value of r/i 23 given by equation ([T])- Firstly, when the tertiary star and the primary star closely 
orbit each other, the tertiary might prevent the primary from formation a PN by causing a too early 
envelope removal. This reduce the denominator in equation ([T]). Secondly, when the higher 
multiple systems are considered the probability of the AGB primary star to interact with a binary 
system rather than with a single secondary star increases. If for example, third of the quadruple, 
quintuplet, and sextuplet stellar systems are taken to be relevant to the present study, then 

|iV(3) + iV(4) + iV(5) + iV(6)]/3 ,,0.26/3 

~/V(2) + Af(3) + Af(4) +/V(5) + Af(6) ~ 0.51 ' ' ' 
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Thirdly, in some systems one of the stars of the tight binary system might be a WD. This is the 
ease when the initially most massive star is in a tight orbit, and it forms a WD without destroying 
the tight binary system. The initially seeond most massive star is in a larger orbital separation 
from the tight binary system. When it evolves on the upper AGB it might swallow the tight binary 
system, now composed of a WD and a main sequence star. 


Two considerations reduce the fraction of PNe shaped by triple-stellar systems relative to the 
value of ? 7 i ,23 given by equations ([1]) and ([21). Firstly, some PNe can result from solar-like stars. 
These have lower mu ltiplicity frequency and companion frequency than stars with initial mass in 
the range 1.5 — 5 Mq (iDuchene & Krausll2013h . Secondly, some PNe might be shaped by planets. I 
assume here that single star do not form bright PNe. But planets might shape PNe as well, adding 
to non-spherical PNe, hence increasing the denominator in equations O and Q. On the other 
hand, one then needs to consider AGB stars that swallow planetary systems. Namely, instead of a 
tight binary system, the AGB swallows a tight planetary system. Such systems can also shape PNe 
simpler to triple-stelar systems if the planet is massive enough. 


Taking these considerations and the large uncertainties I estimate that about one in ten to one 
in six 10% — 17%) of PN progenitors that evolved through a binary interaction, actually went 
through an interaction with a tight binary system. Crudely, about one in eight PNe shaped by 
binary systems was actually shaped by a tertiary system. 


In this work I consider mainly cases where the tight binary system does not survive the CEE. It 
is not crucial for the general results, but it is my expectation that in most, but not all, cases of CEE 
only two stars survive as a bound binary system. One of the two stars in the tight binary system, 
mostly the tertiary star, is expected to merge either with the secondary star or the primary core, or 
be ejected from the system. In a minority of cases all three stars are expected to survive. This 
hypothesis that most tight binary systems do not survive, should be examined by detail numerical 
simulations. 


One case wh ere the tight binary system seems to have survived the CEE is the PN SuWt 2 
(lExter et al.l 120101) . A possible example for evolution where the tight binary system did not sur¬ 
vived the CEE, but all stars survived and stayed bound is the triple pulsar system PSR J0337-I-1715. 
The pulsar is orbite d by two low ma s s whi te dwarfs having periods of 1.63 day and 327 day 
(IRansom et al.ll2014h . ISabach & Sokeri (120ISh proposed the following scenario for this system (for 


an alternative scenario see ITauris & van den Heuvell 120141) . A tight binary system of two main 
sequence stars was tidally and frictionally destroyed inside the envelope of a massive star that 
formed a PN. Only later in the evolution the central massive ONeMg white dwarf experienced an 
accretion-induced collapse and formed the neutron star. That scenario includes a new ingredient 
of a binary system that breaks up inside a CE. As well, the scenario employs an efficient envelope 
removal by jets launched by the main sequence stars immersed in the giant envelope, and the GEE. 
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2.2. Energy considerations 


Two types of energy ean be assoeiated with the tight binary system. The first is the energy 
that is need to breakup the tight binary system, the breakup energy 




_ 1 GM 2 M 3 

23 —- 

2 azAMS,23 


(3) 


The seeond one is the energy the binary system releases if it merges, the merger energy 


^m,23 


1 GM 2 M 2 

2 R 2 


(4) 


where R 2 is the radius of the seeondary star. The amount of energy released depends on the 
response of the seeondary star to the aeeretion of the tertiary star, so the value given above is an 
estimate. 


More eonstruetive is to define dimensionless quantities. Let Of^cE be the final orbital separa¬ 
tion of the see qndary sty from the eore of the prim a ry star aeeording to the energy preseription of 
the CEE (e.g., Webbink 1984 : Tauris & Dewi 2001 : Ivanova et al. 2013 ) that states 


_ 1 GMcoveM2 

■^env — 77 ? 

^ Ctf,CE 


(5) 


where E^env is the binding energy of the giant envelope, ctcE is the CE-a parameter, and Mcove is 
the mass of the giant eore. The two dimensionless ratios are 


t _ -^b,23 

S6 — 7-, 


0.5 


M. 


0.5M„ 


/ 2 af^CE A 

VaZAMS,23/ 


/ Qce \ 
V 0.5 J 


-1 


and 


^ _ -E'm,23 

J-L/ctn-KT 




0.5M,, 


flf.CE 

5i?2 


/ Qce \ 
V 0.5 J 


-1 


( 6 ) 


(7) 


Typieal values of Mcore — O.OMq, Of^cE ~ 2.5Rq (final eore-seeondary orbital period of about 


half a day), R 2 = ORRq {M 2 
the above sealing. 


O. 5 M 0 ), azAMS ,23 ~ 5 i? 0 , and M 3 ^ 0.1 — O. 5 M 0 were used for 


Equation shows that the breakup of the tight binary requires a non-negligible amount of 
energy. This eomes on the expense of the kinetie energy of the seeondary or the tertiary star; one 
star of the tight binary moves to larger orbit, or even eseape the system, and the other star dives 
further in. 


Equation|7]shows that a merger of the tight binary system ean in prineiple ejeet a large fraetion 
of the giant envelope, and at high speeds. The merger oeeurs on several times the dynamieal time 
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of the tight-binary system, whieh is mueh shorter than the orbital period of the triple-system. This 
implies a very asymmetrieal mass ejection to one side. If the two orbital planes, of the tight binary 
system and of the triple system, are parallel, the ejection process maintains the mirror symmetry 
about the orbital plane. But if the two planes are inclined to each other, no axis or plane of 
symmetry are expected in the mass ejection due to the tight-binary merger. A ‘messy’ expanding 
shell is expected to be formed. 

Since most PN central stars have masses in the range Mcore ~ 0.6 — IMq, equations [ 6 ] and |7] 
with the initial constraints M 3 < M 2 < Mzams,i show that on average the role of the tight-binary 
system is expected to increase with primary mass, as the mass of M 3 can be larger. In the case of 
the tight-binary merger process even a tertiary of M 3 0.05 — O.lMcore can play a significant role. 
This includes a brown dwarf tertiary. 


The energy of the merger can be carried by jets, or a wide outflow, launched by the accretion 
disk or belt formed around the secondary star from the destroyed tertiary star. Another source of 
energy, that can exist also in the case of a single secondary star (namely, a binary system), is the 
accretion of giant stellar envelope mass onto the secondary (or tertiary) star. The accretion disk or 
belt can lunch jets t hat facilitate the removal of the giant envelope. The energy carried by the jets 
can be expressed as (ISokerll2014h 


^jets 


-^jets ^ -^acc ®f,CE 1 
Eenv -^core -^2 OiCF, 


Mace A / /«ce\ ^ 

O.lMcore) \^) ’ 


where M^cc is the mass accreted by the secondary (or tertiary) star. 


( 8 ) 


3. OUTCOMES 
3.1. Tight binary in a wide orbit 

3.1.1. Evolution 


The first scenario for a giant star orbited by a tight binary system evolves no CEE or GEE of a 
giant with a tight binary system, but rather accretion by the tight binary system from the dense wind 
of the evolved giant star. It was studied in the past (ISokenbOOdI where more details can be found), 
and it is updated here, as the PNe listed then might not be compatible with new expectations. 


In that study it was found that when the orbital plane of the accreting tight binary system 
and the orbital plane of the triple system (the plane of the tight binary and the giant star) are not 
parallel to each other, the accreted mass onto one or two of the tight binary system components has 
high specific angular momentum. Eor a large fraction of triple-stellar systems accretion disks are 













- 8 - 


expected to form around one or two of the tight binary stars, and luneh jets. The axis of jets will 
be almost parallel to the orbital plane of the triple-stellar system. One jet is blown outward relative 
to the wind, while the other jet passes near the mass-losing star, and is more likely to be deflected 
and slowed down by the giant wind. This flow of asymmetrical jets breaks the mirror symmetry of 
the descendant PN. Namely, the two lobes or two opposites blobs on the two sides of the equatorial 
plane will not be equal. The long-period orbital motion will lead to a departure from axisymmetry, 
i.e., there is no symmetry around the angular momentum axis (the orbital motion itself does not 
break the mirror symmetry about the equatorial plane). 

At the end of the AGB phase when the mass-loss rate is very high, an accretion disk might 
form in tight binary systems up to ~ 400 — 800 AU from the giant star. This leaves a space for 
a fourth star to orbit the AGB star at a close orbit of few AU. The closer stellar companion to 
the AGB star strongly interacts with the AGB star and leads to the formation of an elliptical or a 
bipolar PN. Again, the AGB star itself engulfs only one star, if at all. As we see next, it seems that 
the fourth star is required to form a PN. 


3.1.2. The descendant PNe 


All three stars survive. The triple-system orbit increases due to the mass lost from the AGB 
stars. It can be by a factor of ~ 2 — 3, depending on the initial and final masses. The orbit of the 
tight binary system decreases by a small amount as a result of mass accretion. 


The following PNe were listed in 2004 (ISokeiil2004h as possessing morphologies compatible 
with the expectation from a far-away tight binary system that accretes mass from the wind and 
launches jets. 

He 3-1357 (PN G331.3 -12.1). This PN has a dense ring and small lobes protrude from the nebula 
teobrowsky et al.ll 19981) . However, no single axis of symmetry, nor a point-symmetry, and nor a 
plane of symmetry, can be defined for this nebula. 

PN NGC 6210 (PN G043.1+37.7). This PN, as anoth er example of a PN that was proposed to 
have been shaped by a triple-stellar system (ISokeiil2004t) . is a ‘messy’ PN, with a general elliptical 
structure with unequal sides, blobs, filaments and two pairs of opposite jets protruding from the 


2009: I will return to this PN 


main ‘messy’ (irregular) shell (e.g., IBalicld 119871: IPottasch et al. 
below). __ 

IC 21 49 (PN G166.1+ 10.4) ivazquez et al.ll2002li : Ml -59 (PN G023.9-02.31 iManchado et ^ 


1996 : ISahai et al.ll201 ih : NGC 1514 ( PN G165.5-15.2) (IHajian et al.ll 19971: IMuthu & Anandarao 


2003 ): NGC 6886 (PN G 060.1-07.7) ( Sahai et al. 201 lb . These four PNe were also in that list or 


PNe shaped by triple or quadruple star systems. 
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However, in the years since 2004 a picture has emerged that for a bright P N to be formed 


the mass loss from the AGB star must be enhance b y a strong binary interaction (IDe Marco et al 


20041: ISoker & Subagll2005l:lMoe & De Marcoll2006l) . As such, it is less likely that the above listed 


PNe are compatible with a far-away triple system alone. They might be more compatible with the 
engulfment of a tight binary system (as suggested below for NGC 6210). We therefore need to 
consider one of two scenarios. 


In the first scenario the tight binary system is close to the AGB star, and accretes mass from the 
AGB envelope via a Roche lobe over flow rather than from a wind. The binary system is massive 
enough to maintain the triple system in synchronization and against the Darwin instability, and 
the triple-stellar system survives the entire evolution with no CEE of GEE phases. The outcome 
is a point symmetric bipolar PNe with a triple star system inside: a tight binary system orbiting 
the remnant of the AGB core at ~ 1 AU. In the case that the two orbital planes are parallel, 
the morphology of the descendant PN will be more or less as expected in binary evolution, but 
possibly with a wide and pronounced jets’ precession. When the two orbital planes are inclined to 
each other, a point-symmetric PN with complicated structure that is hard to predict will emerge. 


In the second scenario there is a quadruple-stellar system. The fourth star is closer to the 
AGB star than the far-away tight binary system. Eor the tight binary system to make a mark on 
the descendant PN, it should be closer than 100 AU. I expect the tight binary system to add an 
opposite pair of bullets or two jet-inflated small bubbles. As well, the motion around the center of 
mass of the quadruple systems leads to departure from axisymmetry. Since we are now consider¬ 
ing quadruple systems, such PNe are very rare. It is hard to define the exact characteristics of the 
descendant PNe, and more so to find one. Eet me give one example that might be compatible with 
a quadruple system with a wide tight binary syste m __ 

NGC 6578 (PN GOlO.8-01.8) ( Palen et ^2002 : Sahai et al. 2011 ). This is an elliptical PN con¬ 
taining parts with not well define symmetry dsahai et al. 2011 ). The PN is presented in Pig. [T] It 
has two pairs of opposite unequal structures. One pair has a small lobe on each side of the center, 
but they are unequal in size and shape. Outside there are two opposite bullets, marked with the red 
line on the figure. The other pair has dense structures, one is an arc and the other is a ‘hand fan’ of 
filaments. The two axes are inclined to each other. It is hard to tell at this stage which axes is due 
to the closer fourth star, and which is from the far-away tight binary system. 
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Fig. 1.— An image of NGC 6578 from lSahai et al.l (1201 ih . This morphology might be eompatible 
with a tight binary system in a wide orbit. Aeeording to the proposed seenario, both the eloser 
fourth star and the far-away tight binary systems, eac h launehed a pair o f jets. The two axes of the 
pairs are inclined to each other. One axis marked by ISahai et al.l (1201 lib in dashed-white, and the 
other added here as a red solid line. 


3.2. Ejecting a star 

3.2.1. Evolution 


In the interaction of three point masses one possible outcome is that the lowest mass object is 
ejected from the system, and the two other bodies form a bound binary system. The case where the 
primary star has an envelope of the size of the orbital orbit and the tight binary system breaks up 
inside the giant envelope has not bee n studied numerically. T he breakup of the binary system was 
discussed somewhat sp eculatively by ISabach & Sokeri (120ISh . as summarized in section [2T] here. 
Sabach & Sokeri (12015t) were interested in the case when the two orbital planes are parallel to each 
other. Here we refer to inclined planes as well. 

The initial system is composed of a tight binary system with an orbital separation of azAMS ,23 ~ 
several x lOi?©, that orbit a primary star of mass Mzams,i — 2 — 6Mq, at a triple system orbital 
separation of azAMS,i 23 ~ several x azAMS ,23 ~ 1 — 5 AU. The range of the primary mass is 
dictated by the condition that its maximum radius on the RGB will be much smaller than its max¬ 
imum radius on the AGB when engulfment supposes to take place. T he ratio of the two orbital 
separations ensures dynamical stability (e.g. iMardling & Aarsethll200lh before engulfment occurs, 
and a marginal stability when engulfment starts after the triple system orbit has decreased. 


At the upper AGB the primary star expands. Due to tidal forces the tight binary system spirals- 
in toward the primary star. Now the ratio of the triple system orbital separation to the tight binary 





















- 11 - 


orbital separation decreases, and the system might become dynamical unstable. If the tight binary 
system breaks-up before engulfment occurs, one star, most likely the lowest mass one, is ejected 
from the system, and we enter a regime of binary interaction. 

In the case that the tight binary system does not breakup before engulfment, the two stars start 
to orbit inside the giant envelope. In these systems the orbital separation of the tight binary system 
is large, of the order of magnitude of the density scale hight in the giant envelope. Hence, initially 
only the star closer to the primary is inside the envelope. This star might accrete mass, forms an 
accretion disk, and lunches two opposite jets. At the same time the other star is tens of Rq outside 
the giant photosphere. In the scenario discussed in this subsection, the jets remove the envelope 
outside the binary system. Namely, the system is in the GEE phase. Each of the two stars in the 
tight binary system spends only part of the orbit inside the giant envelope. The jets are expected 
to exist during part of the orbital period of the tight binary system, only when the accreting star(s) 
is(are) inside the envelope. 


Only the star inside the envelope suffers a strong gravitational drag. The other star feels strong 
tidal forces. This situation changes after about a half orbital period, and so on. This differential 
gravitational drag on the two stars and the tidal forces acti ng on the star that is outs ide the envelope, 
might change the outcome in two ways, as speculated bv ISabach & Sokerl (1201511 . Eirstly, instead 
of being ejected from the system, the star that was supposed to escape might end at a large, but 
bound and stable, orbit. Secondly, when the masses ration ^32 = M 3 /M 2 is not much below 1, 
then in some cases the star M 2 > M 3 might end up at a larger orbital separation. 


The break-up is expected to take place when the binary system is in the outer parts of the 
AGB star, in a triple system orbital separation of 0123 > lOOi?©. After the ejection of one star 
from the envelope, the two remaining stars of the inner binary system, either the primary with 
the secondary or the primary with the tertiary, continue in a binary interaction. The inner binary 
system can continue the GEE, hence the two stars survive, or enter a CEE. As in binary interaction 
evolution, the CEE can end with a close binary system, or in a merger. 


One of the four following types of a bound system can be found at the center of the descendant 
PN in cases where the tight binary system breaks apart. (1) A single star. This is the case if one 
of the tight binary stars was ejected and the other merged with the core. (2) A close binary system. 
This is the case if one of the tight binary stars was ejected and the other survived the CEE or the 
GEE phase. (3) A wide and eccentric binary system. This is the case if one of the tight binary stars 
was thrown to a large eccentric orbit, and the other merged with the core. (4) A triple system of a 
wide and eccentric star, and a much closer binary system. This is the case if one of the tight binary 
stars was thrown to a large eccentric orbit, and the other survived the CEE or the GEE phase. If 
the tight binary syster n survives the CE E, the central system will be of a post-AGB star orbited by 
a tight binary system (lExter et al.ll2010h . 
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3.2.2. The descendant PNe 


The descendent nebula possesses a strueture of two non-spherieal struetures resulting from 
two mass loss phases with a relative displaeement between them. The two intensive mass loss 
phases are before and after the breakup of the tight binary system. Before the breakup the mass 
loss oeeurs around the eenter of mass of the triple systems. Jets are formed as the tight binary 
system aeeretes from the AGB envelope via a Roehe lobe overflow. After the ejeetion of one of the 
stars of the tight binary system, the mass loss eontinues around the eenter of mass of the remaining 
binary system. This system moves to the opposite direetion of the ejeeted star. Therefore, the 
nebular part formed after the breakup has a general velocity relative to the outer nebula. This 
descendant PN has a displacement from axisymmetry. If the two orbital planes are not parallel to 
each other, the system will not have even a mirror symmetry about any orbital plane. In cases 
where the tight binary system survives, only the first of these two phases will take place, and the 
descendent PNe will be less co mplicated, althoug h the departure from sphericity can be very large, 
e.g., having an equatorial ring (lExter et al.ll201flt) . 


Let us consider one example out of a very large parameter space. Say at breakup the binary 
system is at 0123 = 1 AU. The masses are Mi = 3Mq, M 2 = O.SMq, and M 3 = 0.3Mq. 
Let us take the tertiary star to become unbound and ejected from the system with a terminal ve¬ 
locity of a fraction x of the escape speed from the system, Vej ,3 ~ 40(x/0.5) km Lrom 
momentum conservation the remaining binary system moves with respect to the center of mass at 
V 12 ~ 3(x/0.5) km s“^. Lor a nebula expanding at ~ 10 — 20 km this leads to an observable 
departure from axisymmetry. 


The PNe are not expected to be too ‘messy’, but might posses no symmetry axis and no plane 
of symmetry. Possible structures are like those of the following PNe. (One should be aware that 
some morphological features can be accounted for also in binary interaction, with no need for 
triple-stellar evolution.) 

IC 2553 (PN G285.4-05.3). It has a point-symmetric structure, but with a displacement of differ¬ 
ent structures (ICorradi et al.ll2000l) . It is possible that a triple-stellar system caused this structure. 
M 1-26 (PN G358.9-00.7). On top of the point-symmetric structure th is PN possesses depart ure 
from axisymmetry along each line that connect two opposite features dSahai & Traugerlll998h . I 
suggest here that a triple-stellar system shaped this PN. An unbound low mass main sequence star 
might exist near the center. 


M 1-6 (PN G211.2-03.5) . This PN has two lobes, and a smaller and de nser inner region dSahai et al. 


201 luOtsuka et al .1120 141) . Its image is presented in Pig. [2] taken from lSahai et al.l d201 ih . This PN 


has neither axisymmetry nor mirror symmetry. Its structure is compatible with the expectation of 
a triple-stellar system shaping discussed in this subsection. 

M 3-1 (PN G242.6-11.6). Has general departures from axisymmetry and from mirror-symmetry 
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PK211-03#1 



L,c,bcr(o),*,h 



Unequal sides caused by ejection of 
a tertiary star. 


Fig. 2.— An image of M 1-6 from ISahai et al.l (1201 ih . I suggest that the departure from both 
mirror-symmetry and from axisymmetry might have been eaused by the ejection of the tertiary star 
from the tight binary system. Despite these deviations, this PN is not as ‘messy’ as expected when 
the tight binary stars merge. 


as those in M 1-6 (ISchwarz et al.lll992h . 


3.3. Drowning star 

3.3.1. Evolution 

In this evolutionary route the triple-stellar system starts more or less as in the stellar ejection 
one, and the tight binary system breaks-up as a result of dynamical instability. However, the dy¬ 
namical instability results in one star moving toward the giant envelope, instead of to the opposite 
direction, and with a low specific angular momentum. The dynamical interaction ‘drowns’ one of 
the stars of the tight binary system in the AGB envelope. This can happen before the tight binary 
system touches the envelope, or after. 

We now have a triple-stellar system composed of a CEE of the AGB star and the drowning 
star of the tight binary system, with a third star orbiting at a distance of Oout ~ 1 AU. There are 
four types of possible outcomes, depending mainly on the relative masses of the three stars. (1) If 
the two stars of the tight binary system are relatively massive, then the drowning star manages to 
eject the AGB envelope. The outcome is a triple-stellar system of a star orbiting the core-remnant 
at Oin ~ 1 — lOi? 0 , and the other one orbiting the central binary system at Oout ~ 1 AU. (2) If the 
drowning star is not massive enough to eject the AGB envelope, it merges with the AGB core. The 
other star can then survive, either with no more interaction with the AGB envelope, or experience a 
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GEE phase. The surviving star ends at a large orbital separation of ~ 1 — 3 AU, depending on how 
mueh mass was ejeeted from the AGB star during the merger process of the core with the drowning 
star. (3) The drowning star does not eject enough AGB envelope mass, and the other star enters a 
CEE phase. This star survives and ends like in a regular binary interaction. (4) The drowning star 
does not eject enough AGB envelope mass, and the other star enters a CEE phase. This star is not 
massive enough, and it merges with the core. The outcome is a single star system. 


3.3.2. The descendant PNe 


The main difference in the morphologies of the descendant PNe in the present evolutionary 
route to those described in sections 13.21 and 13.41 is that no large, or not at all, departure from 
axisymmetry, or more generally no departure from point-symmetry, is expected. Crudely, the 
evolution can be describes as composed of three intensive mass loss episodes, {i) The tight binary 
system accretes mass from the AGB envelope via a Roche lobe overflow. It is expected to launch 
jets and form a bipolar structure, {ii) A CEE of the drowning star inside the AGB envelope. Mass 
loss with concentration toward the equatorial plane might take place, {in) A CEE, or a GEE, or a 
Roche lobe overflow interaction of the other star with the AGB star takes place. Another bipolar 
structure and/or another equatorial mass loss event are expected. If the initial plane of the tight 
binary system and the triple-stellar system are not parallel to each other, the two pairs of lobes 
and/or the two equatorial mass ejections will not be parallel to each other. A large scale point- 
symmetric PN is expected. 


The following PNe have morphologies that might be compatible with the expectation of this 
evolution. 

IPHA S PN-1 This PN conta ins two pairs of lobes, one inside the other, and highly inclined to each 


other (IMampaso et al. 


2006h. There is no large scale departure from point-symmetry 


PN G358.9+03.4 dSahai et al.ll201lh . This PN contains two pairs of lobes highly inclined to each 
other, and one pair is larger than the other. As evident from Pig. |3]no large departure from point- 
symmetry is observed. _ 

Frosty Leo (IRAS 09371+1212) ( Castro-Carrizo et al.12005 1. This is aproto-PN, that like PN G358.9+03.4 
contains two pairs of lobes highly inclined to each other, and one pair is larger than the other. The 
small pair of lobes is actually a complicated structure, where more small lobes are observed. This 

complicated structure might come fr om a triple-stellar evo l ution. __ 

NGC 6302 (PN G349.5+01.0) (e.g., iMeabum et al. 2005; Szvszka et al. 2011 1. It has a well de¬ 
fine bipolar structure, but with another pair of lobes protruding from the main bipolar structure. 

The axis of the protruding pair is inc lined to that of t he bip olar nebula, and the two lobes of the 
protruding pairs are not of equal size. ISoker & Kashil (120121) proposed that the bipolar structure of 
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Fig. 3.— An image of PN G358.9+03.4 from lSahai et al.l (1201 ih . It presents a PN with two pairs 
of lobes highly inelined to eaeh other and no large seale departure from point-symmetry. It might 
be eompatible with the evolutionary route of a drowning star. 


NGC 6302 was formed in an ILOT event. 


3.4. Tight binary merger 

3.4.1. Evolution 

In this ease the two stars enter a CEE or a GEE phase. The gravitational drag on the two stars 
of the tight binary system is more influential than the three-body dynamieal instability, and the two 
stars merge. Aeeretion of mass by one or two of the tight binary stars also deereases the orbital 
separation. Most likely they merge well inside the envelope. As evident from equation ([7]), the 
energy released in the merger proeess is large. A large fraetion of one side of the envelope ean be 
ejeeted at high veloeities, and a binary system survives: the primary and the seeondary with the 
mass it aeereted from the destroyed tertiary star (only a fraetion of the destroyed tertiary star is 
aeereted by the seeondary star). 

Eet us eonsider one example from the large parameter spaee. Say a seeondary of mass and 
radius of M 2 = O.SMq and R 2 = 0.5i?Q, respeetively, aeeretes Mace = O.IMq of the destroyed 
tertiary star. By equation dH) the released energy is ^^,23 ~ 2 x 10^^ erg. 
lasts for about the viseosity time of the aeeretion disk. By equation (2) from 
the viseosity time is about two weeks. If this energy is ehannelled to ejeet, say, Menv-ej = O.SM© 
of the envelope mass, the initial veloeity of the ejeeted envelope mass is Uenv-ej ~ 200 km s“^. 
This is the eseape speed from lOi?© from the eore. So if the merger oeeurs at a radius of 


The aeeretion proeess 


Kashi & Sokerl ( 2015 ) 
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> 2 Oi? 0 , the ejected envelope chunk is ejected at a relatively high speed. 


The fast ejected envelope mass can collide with previously ejected mass, and kinetic energy 
is channelled to thermal energy and radiation. An energy of ~ 10^^ erg can be radiated in sev¬ 
eral weeks. These are close to the radiated energy and time scale of the ILQT of th e AGB star 
NGC 300 OT2008-1 (NGC SOOOTi lMona^boOsilBond et aklbood : [Berger et aklboo^ . So one of 
the outcome of this merger of the tight binary system can be an ILOT event in an AGB star. 


The large released energy facilitates the envelope removal, and the final orbital separation of 
the binary system, the core and secondary star, can be relatively large. As a large mass of the 
envelope is ejected with high momentum to on side, the surviving binary system composed of the 
giant and the secondary star moves to the other. In the example above, let us take the terminal 
velocity of the one-sided ejected mass Menv-ej = O.SMq to be 100 km Let the rest of the 
primary mass be 3.5Mq, which together with the secondary is 4M0. Although the surviving 
system is eight times as a massive as the ejected mass, its velocity relative to the the center of mass 
is lower than 100/8 = 12.5 km s“^. This is because the ejected envelope mass is not a point mass, 
but rather it spreads over a large solid angle. Overall, the surviving binary system might move at 
~ 5 — 10 km relative to the center of mass. The ejection of a large envelope fraction to one 
side, and the motion of the surviving system, that continues to lose mass, to the other side results 
in a ‘messy’ (irregular) PN. 


3.4.2. The descendant PNe 


The descendant nebula has no axisymmetry, and in some cases not even a mirror symmetry 
about the orbital plane. A binary system survives, but it is displaced from the nebular center, 
if a center can be defined at all for the expected ‘messy’ PN. It is expected to have one side 
much brighter than the other. The bright massive side can be larger or smaller than the faint size, 
depending on later evolution, e.g., acceleration by a fast wind. 


The following PNe have morphologies that are compatible with the expectation of a tight 
binary system merger inside the progenitor envelope. 

Sh2-71 (PN G035.9-01). This bipol ^ PN is not really a ‘messy’ PN, but it lacks a well define 
symmetry axis or a symmetry plane (IMocnik et a1.ll2015b . It is compatible with the expectation 
from the scenario discussed here when the energy released in the tight binary merger is not large. 
Abell 46 (PN G055.4+16.0). In this PN the mass distribution departs from any symmetry. The 
companion to the central star has a mass of 0.15Mo (ICorradi et alJl2015b . 

NGC 6210 (PN G043.1+37.7). This was mention bv lsokeij(l2004ll as descendant of a triple-stellar 
system with tight binary being at a wide orbit (section [3. 1.21 above). As evident from Fig. HI this is 
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A messy nebula: 
One-sided mass 
ejection and two 
unequal-iobes pairs 


Fig. 4.— An image of NGC 6210 (from HST site), presenting a ‘messy’ PN that is compatible 
with the expectation from a tight binary system that launched jets and merged. Jet launching can 
occur after the merger process as well. 


indeed a ‘messy’ PN (iBalickil 19871: IPottasch et al.M2009t) . and its morphology is better compatible 
with the expectation from a tight binary merger. The two unequal-lobes pairs might hint that the 
tight binary system launched jets before it merged. 

Hen 2-11 (PN G259.1+00.9). It posses a departure from axisymmetry and mirror symmetry, 
although these deviations are not large, and its cen tral star is a binary system with a K-type main 
sequence companion in a 0.609 d orbital period (I.Tones et al.l 12014t) . The departure from pure 
axisymmetry and mirror symmetry might be compatible with a process of tight binary merger with 
a low mass M-type tertiary star or with a brown dwarf. 

Hen 2-459 (PN G068.3-02.7). It has no lobes. One can define a symmetr y plane, where the mass 
in one side is more extended, and in the other side there is a dense arc dSahai et al.ll2011t) . It is 
presented in Fig. [5l __ 


He 2-71 (PN G296.4-06.9) dSahai & TraugedI 199811: He 2-10 5 (PN G308.6-12.2) (ISchwarz et al. 


1992h : He 2-113 (PN G321.0+03.9) (ISahai et al.ll2000Ll20111) . These are other PNe that lack both 


pure axisymmetry and mirror symmetry, although they have one or more pairs of lobes; some cases 
the two lobes are not equal. 


4. SUMMARY 


I conducted a preliminary study of the possible outcomes from the strong interaction of PN 
progenitors with a close tight binary system. I estimated that about one out of eight PNe that 
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A bright arc to one side couid have been ejected 
by the merger of the tight binary system inside the 
AGB steiiar envelope. 


Fig. 5.— An image of Hen 2-459 from ISahai et al.1 (1201 ih . The pronouneed dipole morphology 


is eompatible with the expectation of mass ejected in a merger process of the tight binary system 


inside the AGB envelope. 


evolved through binary interaction, actually experienced triple-stellar evolution, including brown 
dwarfs tertiary stars (eqs. [I]and[2l). Two types of energy can be associated with the tight binary 
system, the energy that is required to break it up (eq. [3]), and the gravitational energy that is released 
if the two stars of the tight binary system merge (eq. IH). These are not negligible relative to the 
binding energy of the envelope. In particular, the energy released by the merger process can be 
large and leads to ejection of a large portion of the envelope in one direction and possibly with 
high velocities. If part of the kinetic energy is channelled to radiation, an ILOT event of an AGB 
star might take place. 

I also speculated on the possible morphologies of the descendant PNe of triple-stellar evo¬ 
lution. One should be aware that at this stage the affiliation of specific PN morphologies to the 
different triple-stellar evolutionary routes is highly uncertain, and is qualitative and objective. As 
well, many features, like departure from axisymmetry, might be caused also by a binary interaction 
with no need for a triple-stellar system. However, binary interactions are not expected to lead to a 
very large departures from axisymmetry, and are not expected to cause any departure from mirror 
symmetry. 

I concentrated on four evolutionary routes. The first one is of a tight binary system that never 
enters a CEE or a GEE phase, and survives to the PN phase (section ITTI) . To form a bright PN 
either the tight binary system must be close to the AGB star, or a fourth star closer to the AGB 
star is required. In the first case a Roche love over flow is likely to take place, resulting in massive 
jets launched by the tight binary system, that form bipolar PN. In the second case the fourth star 
interacts strongly with the AGB star. Evolutionary routes with a quadruple-stellar system are 
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rare, and it is hard to picture the expected morphology of the descendant PNe of such systems. I 
speculated, non the less, on something like the morphology of NGC 6578 that is presented in Fig. 

m 


In the second evolutionary route the tight-binary system breaks-up as a result of the gravity of 
the primary star (section This occurs when the tight binary system comes closer to the AGB 
star as a result of tidal forces, and a GEE takes place. One star of the tight binary system, more 
likely the tertiary star, is ejected from the system or is thrown to a large and eccentric bound orbit. 
This causes the binary system that is left, of the AGB star and the secondary star, to move in the 
other direction. This leads to departure from axisymmetry of the descendant PN. If the two orbital 
planes, of the tight binary system and of the triple-stellar system, are not parallel to each other, the 
PN will posses no mirror symmetry either. One of the possible expected morphologies might be 
the one of M 1-6 as presented in Fig. [2l 

In the third evolutionary route the dynamical instability breaks-up the tight binary system, and 
throws one star into the AGB envelope with low specific angular momentum; this star ‘drowns’ 
in the AGB envelope (section 13.31) . The evolution can end up with one of four types of stellar 
systems. The descendant PN is expected to have an axisymmetry if the initial equatorial plane of 
the tight binary system is parallel to that of the triple-stellar system, and to posses point-symmetry 
if the planes are inclined to each other. No large scale departure from axisymmetry is expected. I 
speculated that one possible PN morphology is like that shown in Fig. [3l 

In the fourth evolutionary route (section 1X41) the tight binary system merges inside the AGB 
envelope. The large amount of released gravitational energy can eject a large envelope mass in one 
direction. Both before and after merging, and definitely during the merging process, the system 
can launch jets. The PN has no axisymmetry and no point symmetry. If the two orbital planes 
are inclined to each other, the PN will neither posses mirror symmetry. The outcome is a binary 
system with a ‘messy’ PN. Possible PN morphologies are like those shown in Figs. |4]and[51 

In short, the main claim of this paper is that ‘messy’ PNe owe their structure to triple (or 
quadruple) stellar system interaction. 


I did not discuss in details cases where the tight bin ary s ystem survived the GEE and/or the 
GEE phases, a scenario proposed by [Bond et al.l (l2002h and lExter et al.l (120101) . In general the 
morphologies are expected to be similar to some of the cases of the drowning stellar scenario 
(section IXX) . If, for example, the two orbital planes are parallel, then th e outcome can be as in 


binary evolution, but more extreme, e.g., an extreme ring as suggested bv (Bond et al.l (120021) and 
Exter et al.l (120101) . If the two planes are inclined, a point symmetric PNe might be formed, with a 
pair of lobes inclined at an angle different than 90° to a dense equatorial ring or torus. Namely, 
the descendant PNe of surviving tight binary systems are expected to posses axial-symmetry or 
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point-symmetry despite the large departure from spherieal-symmetry. They will not be ‘messy’. 

I benefitted from discussions with Amit Kashi, Erez Michaely, and Efrat Sabach on several 
processes mentioned in the paper. I thank Howard Bond, David Jones, Raghvendra Sahai, and Hans 
Van Winckel, and an anonymous referee for valuable comments and suggestions. I am supported 
by the Charles Wolfson Academic Chair. The Planetary Nebula Image Catalogue (PNIC) compiled 
by Bruce Balick was an essential tool in this study 
(http ://w w w. astro .washington.edu/users/balick/PNIC/). 


REFERENCES 

Akras, S., Boumis, P, Meaburn, J., Alikakos, J., Eopez, J. A., Goncalves, D. R. 2015, MNRAS, 
452, 2911 

Aller, A., Miranda, E. E, Olgum, E., Vazquez, R., Guillen, P. E, Oreiro, R., Ulla, A., & Solano, E. 
2015a, MNRAS, 446, 317 

Aller, A., Montesinos, B., Miranda, E. E, Solano, E., & Ulla, A. 2015b, MNRAS, 448, 2822 
Balick, B. 1987, AJ, 94, 671 

Berger, E., Soderberg, A. M., Chevalier, R. A., et al. 2009, ApJ, 699, 1850 
Bobrowsky, M., Sahu, K. C., Parthasarathy, M., & Garcia-Eario, P. 1998, Nature, 392, 469 
Boffin, H. 2015, 19th European Workshop on White Dwarfs, 493, 527 

Bond, H. E. 2000, Asymmetrical Planetary Nebulae II: Prom Origins to Microstructures, 199, 115 

Bond, H. E., Bedin, E. R., Bonanos, A. Z., Humphreys, R. M., Monard, E. A. G. B., Prieto, J. E., 
& Walter, P. M. 2009, ApJ, 695, E154 

Bond, H. E., & Eivio, M. 1990, ApJ, 355, 568 

Bond, H. E., O’Brien, M. S., Sion, E. M., Mullan, D. J., Exter, K., Pollacco, D. E., & Webbink, 
R. P. 2002, Exotic Stars as Challenges to Evolution, 279, 239 

Bujarrabal, V, Alcolea, J., Van Winckel, H., Santander-Garcia, M., & Castro-Carrizo, A. 2013, 
A&A, 557, A104 

Castro-Carrizo, A., Bujarrabal, V, Sanchez Contreras, C., Sahai, R., & Alcolea, J. 2005, A&A, 
431,979 


- 21 - 


Corradi, R. L. M., Garcfa-Rojas, J., Jones, D., & Rodriguez-Gil, R 2015, ApJ, 803, 99 

Corradi, R. L. M., Gon 9 alves, D. R., Villaver, E., Mampaso, A., & Perinotto, M. 2000, ApJ, 542, 
861 

Decin, L., Richards, A. M. S., Neufeld, D., Steffen, W., Melnick, G., & Lombaert, R. 2015, A&A, 
574, A5 

De Marco, O. 2015, in “Physics of Evolved Stars 2015 - A conference dedicated to the memory of 
Olivier Chesneau”, in press 

De Marco, O., Bond, H. E., Harmer, D., & Eleming, A. J. 2004, ApJ, 602, E93 

De Marco, O., Eong, J., Jacoby, G. H., Hillwig, T., Kronberger, M., Howell, S. B., Reindl, N., 
Margheim, S. 2015, MNRAS, 448, 3587 

De Marco, O., Passy, J.-C., Moe, M., Herwig, E, Mac Eow, M.-M., & Paxton, B. 2011, MNRAS, 
411,2277 

De Marco, O., & Soker, N. 2011, RASP, 123, 402 

Douchin, D., De Marco, O., Erew, D. J., Jacoby, G. H., Jasniewicz, G., Eitzgerald, M., Passy, J-C., 
Harmer, D., Hillwig, T., & Moe, M. 2015, MNRAS, 448, 3132 

Duchene, G., & Kraus, A. 2013, ARA&A, 51, 269 

Exter, K., Bond, H. E., Stassun, K. G., Smalley, B., Maxted, P E. E., & Pollacco, D. E. 2010, AJ, 
140, 1414 

Eang, X., Guerrero, M. A., Miranda, E. E, Riera, A., Velazquez, P E, Raga, A. C. 2015, MNRAS, 
452, 2445 

Garcia-Segura, G., Villaver, E., Eanger, N., Yoon, S.-C., & Manchado, A. 2014, ApJ, 783, 74 

Gorlova, N., Van Winckel, H., Gielen, C., et al. 2012, A&A, 542, A27 

Gorlova, N., Van Winckel, H., Ikonnikova, N. P, Burlak, M. A., Komissarova, G. V, Jorissen, A., 
Gielen, C., Debosscher, J., & Degroote, P 2015, MNRAS, 451, 2462 

Hajian, A. R., Prank, A., Balick, B., & Terzian, Y. 1997, ApJ, 477, 226 

Hillwig, T. C., Erew, D. J., Eouie, M., De Marco, O., Bond, H. E., Jones, D., Schaub, S. C. 2015, 
AJ, 150, 30 

Ivanova, N., Justham, S., Chen, X., et al. 2013, A&A Rev., 21, 59 



- 22 - 


Jones, D. 2015, in “Physics of Evolved Stars 2015 - A conference dedicated to the memory of 
Olivier Chesneau”, in press (arXiv: 1507.05447) 

Jones, D., Boffin, H. M. J., Miszalski, B., Wesson, R., Corradi, R. L. M., & Tyndall, A. A. 2014, 
A&A, 562, A89 

Jones, D., Boffin, H. M. J., Rodriguez-Gil, R, Wesson, R., Corradi, R. L. M., Miszalski, B., & 
Mohamed, S. 2015, A&A, 580, A19 

Kashi, A., & Soker, N. 2015, arXiv: 1508.00004 

Lagadec, E., Verhoelst, T., Mekamia, D., et al. 2011, MNRAS, 417, 32 

Eombardi, J. C., Jr., Proulx, Z. E, Dooley, K. E., Theriault, E. M., Ivanova, N., & Rasio, E. A. 
2006, ApJ, 640, 441 

Mampaso, A., Corradi, R. E. M., Viironen, K., et al. 2006, A&A, 458, 203 

Manchado, A., Guerrero, M. A., Stanghellini, E., & Serra-Ricart, M. 1996, The lAC morphological 
catalog of northern Galactic planetary nebulae. Publisher: Ea Eaguna, Spain: Institute de 
Astrofisica de Canarias (lAC), 1996, Eoreword by Stuart R. Pottasch, ISBN: 8492180609, 

Manick, R., Miszalski, B., & McBride, V. 2015, MNRAS, 448, 1789 

Mardling, R. A., & Aarseth, S. J. 2001, MNRAS, 321, 398 

Martinez Gonzalez, M. J., Asensio Ramos, A., Manso Sainz, R., Corradi, R. E. M., & Eeone, E. 
2015, A&A, 574, A16 

Meabum, J., Eopez, J. A., Steffen, W., Graham, M. E, & Holloway, A. J. 2005, AJ, 130, 2303 

Miszalski, B., Manick, R., & McBride, V. 2015, in “Physics of Evolved Stars 2015 - A conference 
dedicated to the memory of Olivier Chesneau”, in press (arXiv: 1507.07707) 

Mocnik, T., Eloyd, M., Pollacco, D., & Street, R. A. 2015, MNRAS, 451, 870 

Moe, M., & De Marco, O. 2006, ApJ, 650, 916 

Monard, E. A. G. 2008, lAU Circ., 8946, 1 

Montez, R., Jr., Kastner, J. H., Balick, B., et al. 2015, ApJ, 800, 8 

Muthu, C., & Anandarao, B. G. 2003, AJ, 126, 2963 

Nordhaus, J., & Blackman, E. G. 2006, MNRAS, 370, 2004 


-23- 


Otsuka, M., Kemper, R, Cami, J., Peeters, E., & Bemard-Salas, J. 2014, MNRAS, 437, 2577 

Palen, S., Balick, B., Hajian, A. R., Terzian, Y., Bond, H. E., & Panagia, N. 2002, AJ, 123, 2666 

Passy, J.-C., De Mareo, O., Eryer, C. E., et al. 2012, ApJ, 744, 52 

Passy, J.-C., De Mareo, O., Eryer, C. E., et al. 2012, ApJ, 744, 52 

Pottaseh, S. R., Bernard-Salas, J., & Roellig, T. E. 2009, A&A, 499, 249 

Ransom, S. M., Stairs, I. H., Archibald, A. M., et al. 2014, Nature, 505, 520 

Ricker, P. M., & Taam, R. E. 2012, ApJ, 746, 74 

Sabach, E., & Soker, N. 2015, MNRAS, 450, 1716 

Sahai, R., Morris, M. R., & Villar, G. G. 2011, AJ, 141, 134 

Sahai, R., Nyman, E.-A., & Wootten, A. 2000, ApJ, 543, 880 

Sahai, R., & Trauger, J. T. 1998, AJ, 116, 1357 

Sandquist, E. E., Taam, R. E., Chen, X., Bodenheimer, R, & Burkert, A. 1998, ApJ, 500, 909 

Schwarz, H. E., Corradi, R. E. M., & Melnick, J. 1992, A&AS, 96, 23 

Soker, N. 1994, MNRAS, 270, 774 

Soker, N. 1996, MNRAS, 283, 1405 

Soker, N. 1998, ApJ, 496, 833 

Soker, N. 2002, MNRAS, 330, 481 

Soker, N. 2004, MNRAS, 350, 1366 

Soker, N. 2013, New A, 18, 18 

Soker, N. 2014, arXiv: 1404.5234 

Soker, N. 2015, ApJ, 800, 114 

Soker, N., & Harpaz, A. 1992, PASP, 104, 923 

Soker, N., & Kashi, A. 2012, ApJ, 746, 100 

Soker, N., & Subag, E. 2005, AJ, 130, 2717 


-24- 


Soker, N., Zucker, D. B., & Balick, B. 1992, AJ, 104, 2151 
Szyszka, C., Zijlstra, A. A., & Walsh, J. R. 2011, MNRAS, 416, 715 
Tauris, T. M., & Dewi, J. D. M. 2001, A&A, 369, 170 
Tauris, T. M., & van den Heuvel, E. R J. 2014, ApJ, 781, L13 
Tokovinin, A., Thomas, S., Sterzik, M., & Udry, S. 2006, A&A, 450, 681 
Tsebrenko, D., & Soker, N. 2015, MNRAS, 447, 2568 

Van Winckel, H., Jorissen, A., Exter, K., Raskin, G., Prins, S., Perez Padilla, J., Merges, E, 
Pessemier, W. 2014, A&A, 563, ElO 

Vazquez, R., Miranda, E. E, Torrelles, J. M., Olguin, E., Benitez, G., Rodriguez, E. E, Eopez, 
J. A. 2002, ApJ, 576, 860 

Webbink, R. E. 1984, ApJ, 277, 355 

Zijlstra, A. 2015, Revista Mexieana de Astronomia y Astroffsiea, aeeepted for publieation 
(arXiv:1506.05508) 


This preprint was prepared with the AAS lATj^ macros v5.2. 



